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ABSTRACT 

We analyse the kinematics of the entire spectroscopic sample of 99 recently discovered 
high proper-motion white dwarfs by Oppenheimer et al. using a maximum-likelihood 
analysis, and discuss the claim that the high-velocity white dwarfs are members of a 
halo population with a local density at least ten times greater than traditionally as- 
sumed. We argue that the observations, as reported, are consistent with the presence 
of an almost undetected thin disc plus a thick disc, with densities as conventionally as- 
sumed. In addition, there is a kinematically distinct, flattened, halo population at the 
more than 99% confidence level. Surprisingly, the thick disc and halo populations are 
indistinguishable in terms of luminosity, color and apparent age (1-10 Gyr). Adopt- 
ing a bimodal, Schwarzschild model for the local velocity ellipsoid, with the ratios 
a\j:ay:ayj~l:2/3:l/2, we infer radial velocity dispersions of eru=62t 10 kins" 1 and 

150^40 kms -1 (90% C.L.) for the local thick disc and halo populations, respectively. 
The thick disc result agrees with the empirical relation between asymmetric drift and 
radial velocity dispersion, inferred from local stellar populations. The local thick-disc 
plus halo density of white dwarfs is 7i* d +^=(1.9 ± 0.5) xl0~ 3 pc~ 3 (90% C.L.), of 
which Tig WD =l.llg'7xlO~ 4 pc~ 3 (90% C.L.) belongs to the halo, a density about five 
times higher than previously thought. Adopting a mean white-dwarf mass of 0.6 Mq, 
the latter amounts to CSIq^xIO -2 (90% C.L.) of the nominal local halo density. 
Assuming a simple spherical logarithmic potential for the Galaxy, we infer from our 
most-likely model an oblate halo white-dwarf density profile with n(r) oc r~ a and 
a « 3.0. The halo white dwarfs contributes ~ 2.6 x 10 9 M , i.e. a mass fraction of 
^0.004, to the total mass inside 50 kpc (JIwd ~ 10~ 4 ). The halo white dwarf popu- 
lation has a microlensing optical depth towards the LMC of t^, d m 1.3 x 10~ 9 . The 
thick-disc white dwarf population gives r^ D w 4 x 10~ 9 . The integrated Galactic 
optical depth from both populations is 1-2 orders of magnitude below the inferred 
microlensing optical depth toward the LMC. If a similar white-dwarf population is 
present around the LMC, then self-lensing can not be excluced as explanation of the 
MACHO observations. We propose a mechanism that could preferentially eject disc 
white dwarfs into the halo with the required speeds of ^200 kms -1 , through the or- 
bital instability of evolving triple star systems. Prospects for measuring the density 
and velocity distribution of the halo population more accurately using the Hubble 
Space Telescope Advanced Camera for Surveys (ACS) appear to be good. 
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1 INTRODUCTION 

As the most common stellar remnants, white dwarfs (WD) 
provide an invaluable tracer of the early evolution of the 
Galaxy. Their current density and distribution reflects the 
disposition of their progenitor main sequence stars and their 
colors indicate their ages, which are consistent with many of 
them having been formed when the Galaxy was quite young 
(e.g. Wood 1992). 



Until recently, it has been assumed that halo WDs con- 
tribute a negligible fraction of the total mass of the Galaxy 
(e.g. Gould, Flynn & Bahcall 1998). This view is supported 
by the constraint that the formation of ~O.6-M0 WDs is 
accompanied by the release of several solar masses of gas, 
heavily enriched in CNO elements (e.g. Chariot & Silk 1995; 
Gibson & Mould 1997; Canal, Isern & Ruiz-Lapuente 1997). 
Yet local stars and interstellar gas, which only comprise 
~ 10% of the total Galaxy mass, contain only 2-3 percent 
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Figure 1. Color-color diagram of the 
sample of WDs. Shown are the high- 
velocity (>94 kms -1 ; squares) and low- 
velocity WDs (<94 kms -1 ; triangles). 
The velocity cuts in U and V (see 
Sect. 2) are taken as in Oppenheimer et 
al. (2001), for reasons of consistency, 
but assume v r = kms -1 . The cool- 
ing curves are those from Chabrier et 
al. (2000). The upper to lower curves rep- 
resent 0.5, 0.6 and 0.8-M Q WDs with 
pure hydrogen atmospheres. Ages in Gyr 
are indicated along the 0.5-Mq curve. 



of these elements by mass. Based on these metalicity ar- 
guments, most recently Fields, Freese & Graff (2000) have 
argued that £1wd<3 x 10~ 4 from C and N element abun- 
dances, adopting Ho=70kms -1 Mpc -1 . WDs can therefore 
not comprise more than 1-10% percent of the dark matter 
in the Galactic halo within 50 kpc, if f2h a io(< 50 kpc) ~ 0.03 
(e.g. Bahcall et al. 2000) . However, this argument can be cir- 
cumvented if significant metal outflow from the Galaxy, due 
to supernovae, accompanies the formation of these WDs, 
thereby removing most of the produced metals from the 
Galaxy (e.g. Fields, Mathews & Schramm 1997). 

These considerations are important because the MA- 
CHO collaboration (e.g. Alcock et al. 2000) report a fre- 
quency of Large Magellanic Cloud (LMC) microlensing 
events that is suggestive of 8-50% (95% C.L) of the halo 
mass in ~ O.5-M0 compact objects inside a 50 kpc Galacto- 
centric radius. Faint stars and brown dwarfs seem to com- 
prise only a few percent of the Galaxy mass (e.g. Bahcall 
et al. 1994; Graff & Freese 1996a, 1996b; Mera, Chabrier & 
Schaeffer 1996; Flynn, Gould & Bahcall 1996; Freese, Fields 
& Graff 1998) and are therefore unlikely to be responsible for 
these microlensing events. The only remaining candidates 
with similar masses, known to us, are WDs. 

Alternative explanations for these microlensing events, 
based on stellar self-lensing in the LMC have been put for- 
ward (Wu 1994; Sahu 1994; Salati et al. 1999; Aubourg et 
al. 1999; Evans & Kerins 2000). A comparable microlensing 
survey done by the EROS collaboration only claims an up- 
per limit of 40% (95% C.L.) of the halo mass in compact 
objects with mass less than 1 Mq (Lasserre et al. 2000; see 



Afonso et al. 1999 for SMC results). Because the EROS sur- 
vey targets the outer regions of the LMC, in contrast to the 
MACHO survey, this could indeed suggest some self-lensing 
contribution. 

There have also been claims of a number of high proper- 
motion objects in the HDF, which were suggested to be old 
WDs in the Galactic halo at a distance of ~2 kpc (Ibata 
et al. 1999), with blue colors that were consistent with re- 
vised white-dwarf atmosphere models (Hansen 1999). More 
locally, Ibata et al. (2000) found two high proper-motion 
WDs. If the HDF results had been confirmed, they alone 
would have been sufficient to explain the observed microlens- 
ing optical depth towards the LMC. However, recent third- 
epoch observations have shown that these "objects" were 
misidentifications and not moving objects (Richer 2001). 

The topic has become interesting again, though, with 
the recent discovery by Oppenheimer et al (2001) of 38 
nearby, old WDs, with de-projected horizontal velocities in 
excess of 94 km s _1 (see Sect. 2). They concluded that at least 
two percent of the halo is composed of these cool WDs. This 
conclusion has been challenged, however, by Reid, Sahu & 
Hawley (2001), who claim that these WDs could be the high- 
velocity tail of the thick disk. In addition, Hansen (2001) has 
argued that these WDs follow the same color-magnitude re- 
lation as those of the disc, which might be unexpected if 
they are a halo population. 

In this paper, we take the observations of Oppenheimer 
et al. (2001) at face value and examine their conclusions in 
more detail. We assume that the high proper motion objects 
have correctly been identified as WDs and that their inferred 
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distances are statistically unbiased. However, it is important 
that additional observations be performed to validate this. 
We also discuss how to set more precise bounds on the halo 
WD density. 



2 THE WHITE-DWARF SAMPLE 

The high proper motion WDs were selected from digitized 
photographic plates from the SuperCOSMOS Sky Surveyrj 
(Hambly et al. 2001), covering an effectively searched area of 
4165 square degree around the South Galactic Pole (SGP), 
upto b~— 45°. Objects with proper motions between /i=0.33 
and 10.0 arcsec yr _1 were selected, that were brighter 
than R59F=19.8 mag and seen at three epochs in each 
field. From these objects a sample of 126 potential halo 
WDs was drawn based on their reduced proper motion 
(Hn=R59F+5log fi+5) and their color Bj — R. After spec- 
troscopic observations of 69 of the 92 WDs without spectra 
(34 had published spectra), a sample of 99 WDs remained. 
For a full discusion of the sample, its selection criteria and 
completeness, we refer to Oppenheimer et al. (2001). 

Distances to each of these WDs were ascribed using 
the photometric parallax relation derived by Bergeron et al. 
(1997). This was checked in two cases using measured dis- 
tances and the relative uncertainty was estimated at 20 per- 
cent. Using these distances, Oppenheimer et al. (2001) sub- 
sequently performed a velocity selection on the WD sam- 
ple. The proper motions were converted into a velocity in 
Galactic coordinates relative to the local standard of rest, 
(LSR), with U directed towards the Galactic center (GC), 
V along the direction of the circular velocity and W to- 
wards the North Galactic Pole (NGP), see Fig.2. The ab- 
sence of a measured radial velocity (v r ) was compensated 
by assuming that 1^=0 and deprojecting onto the U — V 
plane, ft was then argued that 95 percentP] of the thick-disc 
WDs have velocities in the plane within 94 km s _1 (« 2a t d) 
from the asymmetric drift point ([/, V) = (0, — 35) km s _1 
and consequently all such WDs were eliminated, leaving a 
sample of 38 high-velocity WDs of which 26 are new discov- 
eries and 14 exhibit hydrogen lines and are believed to be 
younger. These WDs have velocities in the Galactic plane 
that are enclosed within a circle of radius 2<Th centered on 
(U,V) = (0,-220) km s _1 corresponding to a non-rotating 
halo distribution. The density of WDs in the magnitude- 
limited sample was then estimated to be ~ 10 -4 Mq pc -3 , 
roughly ten times the expected density (e.g. Gould, Flynn 
& Bahcall 1998) and equivalent to 2 percent of the nominal 
local halo mass density (e.g. Gates, Gyuk & Turner 1995). 

This procedure has been criticised by Reid, Sahu & 
Hawley (2001) who de-project the velocity vector by setting 
the radial velocity (v r ) to zero, thereby reducing the number 
of high velocity WDs. This procedure, however, also places 
several lower velocity WDs outside the 94 km s _1 cut. Reid 
et al. (2001) also note that the high velocity WDs are not 
concentrated around (U, V) — (0, —220) km s _1 as might be 
expected from a halo population and argue that the WDs 
are mostly associated with the high velocity tail of the thick 
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Figure 2. The Galactic coordinate system. The U, V and W 
velocity components point toward the Galactic center (GC), in 
the direction of Galactic rotation (anti-clockwise as seen from 
the North Galactic Pole, NGP) and out of the Galactic Plane, re- 
spectively. The velocity components ti£ and Vi are the projected 
velocity components on the sky, in the Galactic coordinates sys- 
tem (I, b), of a WD with space velocity -uwDi and v r is the radial 
velocity. The vector vq indicates the Solar motion with respect 
to the Local Standard of Rest (LSR). 



disc plus a few WDs from the traditional halo stellar pop- 
ulation. Essentially the same point has been made by Graff 
(2001). fn addition, Gibson & Flynn (2001) identified a num- 
ber of errors in the original table of white-dwarf properties 
and argued that the Oppenheimer et al. (2001) densities are 
further overestimated by a factor 3-10, bringing them more 
into line with the results from other surveys. 

In this paper, we analyse the full sample of 99 WDs, 
not only those 38 with the highest space velocities, thereby 
avoiding potential pitfalls associated with arbitrary cuts in 
the sample and analysing only subsamples. 



2.1 White-Dwarf Distances and Ages 

The photometric parallax method to obtain distances hides 
a potential problem for the oldest WDs (<T0 Gyr). Al- 
though these white dwarfs continue to cool and fade away, 
their Bj — R59F color remains nearly constant at sal. 9 for 
M Bj £l8 (Fig.l; e.g. Chabrier et al. 2000). Blindly applying 
the linear color-magnitude relation given in Oppenheimer et 
al. (2001) to the oldest WDs would therefore lead to a se- 
vere overestimate of their intrinsic luminosity and distance. 
Even if old WDs are present in the sample, it might not be 
surprising that they don't show up in the color-magnitude 
relation by Hansen (2001). WDs that move onto the hori- 
zontal branch of the cooling curve have increasingly overesti- 
mated luminositities and consequently move up in the color- 
magnitude diagram much closer to the sequence of younger 
(<J10 Gyr) WDs. Fortunately, only 5% of the sample of WDs 



have colors Bj — R59F^,1.8 (Fig.l) and it should therefore 
not create a major problem. 

However, it is potentially dangerous to create color- 
magnitude diagrams, using distances that have been derived 
by forcing the WDs to follow the color-magnitude relation 
of relatively young WDs, and from those to conclude that 
the WD are young. This is a circular argument. A better ar- 
gument is to plot the age sequence on a two-color diagram 
(Fig.l; see also Oppenheimer et al. 2001). Again, there is 
no evidence for WDs older than ~ 10 Gyr and the high- 
velocity WDs appear to be drawn from the same popula- 
tion as the low-velocity WDs which are presumably relics of 
ongoing star formation in the disc. This appears to argue 
against there being a substantial white-dwarf halo popula- 
tion (e.g. Hansen 2001), although if there exists a mechanism 
that preferentially gives WDs a high space-velocity during 
its formation/evolution one might not expect a correlation 
between age and velocity. The latter would only be expected 
if velocities increased with time due to slower scattering pro- 
cesses. We return to this in Sect. 7. 

We should emphasize, however, that the implied ages 
and distances are only as good as the assumed underlying 
white-dwarf color-magnitude relation, in this case that from 
Bergeron et al. (1997), and that independent measurements 
of their distances are necessary, preferentially through par- 
allax measurements, even though that will be difficult given 
their large proper motions. 



3 COORDINATES AND KINEMATICS 

For each WD, we transform its equatorial coordinates (a, 5) 

in to Galactic coordinates (Z,i>)Fl Using the distances as- 
cribed by Oppenheimer et al. 2001, we then transform the 
observed proper motion vector in Galactic coordinates into 
a vector (ve,Vb), where ve is velocity component along the 
/-direction and Vb is the velocity component along the the 
fa-direction (Fig. 2). The resulting vector is orthogonal to the 
line-of-sight, i.e. the direction of the radial velocity (v r ), 
which is unknown to us. The transformation from the ob- 
served velocities to Galactic velocities (U,V ,W) with respect 
to the Local Standard of Rest (LSR), is 

U — Uq + v r cos b cos £ — Vb sin b cos £ — ve sin £ 
V — Vq + v r cos b sin £ — Vb sin b sin £ + vi cos £ 
W — Wq + v r sin b + v b cos b (1) 

where (Uq,Vq,Wq) = (+10, +5.2, +7.2) km s" 1 is the Solar 
motion with respect to the LSR (Dehen & Binney 1998b). 
We can use Eq.rtTJ) to transform the velocity distribution 
function f(U,V,W) (§4.6) in to f(vi,Vb,v r ), or visa versa, 
at the position of any WD with a measured proper-motion 
vector and distance. 

Fig. 3 shows velocity vector (V(,Vb) for each WD at its 
corresponding Galactic coordinates. The velocity vectors 
point preferentially in the direction opposite to Galactic ro- 
tation, which reflects the asymmetric drift (§4.2), whereby 
a stellar population that is significantly pressure supported, 

•f We use the transformation in §2.1.2 from BM98. Note the ty- 
pographic error in Zcp=122.932° for the Galactic longitude of the 
Celestial Pole. 



rotates slower than its LSR. If many of the high proper mo- 
tion objects had been misidentied on the photometric plates, 
one does not expect a collective movement of these objects 
in the same direction. 



4 A PHASE-SPACE DENSITY MODEL 

Before introducing a model for the phase-space density of 
WDs in the Solar neighborhood in §4.6, we first (i) describe 
a simple global potential model that we will employ to derive 
some of the characteristic orbital parameters of the sample 
of WDs, (ii) reprise epicyclic theory, in order to clarify the 
discussion and some of the concepts that we will use in this 
paper, (iii) derive some characteristic quantities of the white- 
dwarfs data set, (iv) examine the properties of the white- 
dwarf orbits and (v) discuss some general properties of the 
Galactic disc components. 



4.1 A Toy Model of the Galactic Potential 

For definiteness we assume a spherically symmetric potential 
for the Galaxy 



$(R,z) = (v 2 c /2)-\n(R 2 + z 2 ), 



(2) 



in cylindrical coordinates, where R is the distance from the 
rotation axis of the Galaxy and z is the height above the 
plane. The simplification of spherical symmetry is not incon- 
sistent with observations of the Galaxy and external galax- 
ies, even though a wide range of values for the flattening 
of the Galactic halo is allowed, often strongly dependent on 
the method used to derive it (see Oiling & Merrifield 2000). 
This potential, although a clear simplification, it is to first 
order also consistent with the relatively flat Galactic rota- 
tion curve (e.g. Dehnen & Binney 1998a) and allows us to 
gain insight in the kinematics and properties of the white- 
dwarf sample. The circular velocity v c =d&/d\nR in this 
potential is constant in the plane of the Galaxy (z=0). 

One can subsequently introduce the effective potential 



$ cS (R,z) = <S>(R,z) + ^, 



(■i) 



where L z = Rv,/, is the angular momentum of a WD and 
v,f, is its tangential velocity in the plane of the Galaxy. If a 
particle is "at rest" in the minimum of the effective poten- 
tial, i.e. d<& a fs/dR=0, one readily finds that v^,=v c , hence 
the WD moves on a circular orbit with a constant velocity 
v c =Q(R)R, where Q(R) is the angular velocity. 



4.2 Epicyclic Orbits 

Most stars and stellar remnants, however, do not move on 
perfect circular orbits, but have an additional random ve- 
locity and move on orbits around the minimum in their ef- 
fective potential. If these excursions have relatively small 
anrplitudes, the Taylor expansion of the effective potential 
around this minimum can be terminated at the second- 
order derivative and the stars or WDs will exhibit sinusoidal 
(i.e. epicyclic) motions around the circular orbit. The fre- 
quency of these epicyclic motions are k 2 = d 2 Q e s/dR 2 = 
(R~ 3 d(R 4 Q 2 ) J dR) r s in the tangential and radial directions 
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Figure 3. The white-dwarf velocity vec- 
tors (v£,vi,), with respect to an observer 
near the Sun, projected on the sky in 
Galactic coordinates. Galactic latitude 
increases from —90° to —30° in steps of 
15° (dashed circles). Galactic longitude 
increase anti-clockwise with / = 0° to- 
wards to Galactic center (GC). The ar- 
rows indicates the velocity and direction 
of the WDs. The asymmetric drift of the 
higher velocity WDs (>150 kms -1 ; thick 
arrows) is clearly visible, pointing in the 
direction opposite to Galactic rotation 
(£ = 90°). 



and v 2 = d 2 & c g/dz 2 in the direction out of the plane. The 
radius of the guiding centre is 7? g (see §4.2). 

Because these frequencies depend on local higher-order 
derivatives of the potential, we can not use our simple po- 
tential model, Eqn.(fe|), but have to rely on observational 
constraints (§4.2). One can introduce the Oort's constants 
A and B, whose values have observationally been deter- 
mined to be j4=14.5kms -1 kpc -1 and B——12kms^ 1 kpc -1 
in the Solar neighborhoods (BT87). From this is follows that 
«g = -AB(A-B) = -ABQo, i.e. k = 36 kms" 1 kpc -1 , and 






= 2 



-B 



(A-B) 



1.3. 



(4) 



The frequency of the epicyclic motion is therefore constant 
at 1.3 times its orbital frequency, irrespective of the ampli- 
tude of the orbit around the minimum of its effective po- 
tential. In comparison, the potential given by Eq.(Q) gives 
K /n = a/2 « 1.4. 

Under the epicyclic approximation, stars undergo ret- 
rograde, elliptical, orbits in the plane about their guiding 
centers, which move on circular orbits. The ratio of the ra- 
dial velocity amplitude to the azimuthal velocity amplitude 
is /t/2fi w 0.7, where Q = 27 km s _1 kpc" 1 is the circular an- 
gular frequency (e.g. BT87) and ko is given above. The indi- 
vidual orbits are elongated azimuthally. However, at a given 



8 All subscripts indicate values in the Solar neighborhood 
at Rq. 



point, retrograde motion by stars with interior guiding cen- 
ters that orbit faster than the LSR, must be combined with 
prograde motion executed by stars with exterior guiding cen- 
ters that orbit more slowly. The net result is that the ratio 
of the radial velocity dispersion to the azimuthal velocity 
dispersion is predicted to be eru/av=2f2/««1.5. This pre- 
diction is supported by observations of different stellar pop- 
ulations (e.g. Dehnen & Binney 1998b; Dehen 1998; Chiba 
& Beers 2000). 

Superimposed on this horizontal motion is a vertical os- 
cillation with angular frequency, v — 3.6 H, which is uncou- 
pled to the horizontal motion at this level of approximation. 
The amplitude of this oscillation is determined by the past 
history including scattering off molecular clouds and spiral 
arms, and is measured to be approximately aw/cru~0-5 (e.g. 
Dehnen & Binney 1998b; Dehen 1998; Chiba & Beers 2000). 
The velocity ellipsoid appears to maintain these dispersion 
ratios even for velocities where the epicyclic approximation 
ought to be quite inaccurate. 

At the next level of approximation, we must include the 
radial density variation and this will contribute a mean tan- 
gential velocity, relative to the LSR, at a point, as there are 
more stars with interior guiding centers moving backwards 
than exterior stars moving forwards. This is the asymmetric 
drift, v a —— <V>=v c — Vj,, which is the velocity lag of a 
population of stars with respect to the LSR. This can be es- 
timated from the fluid (Jeans') equations by observing that 
the circular velocity of a group of stars will be reduced by 
a quantity proportional to the pressure gradient per unit 
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Figure 4. The projection of the 99 
white-dwarf velocity vectors on the U— V 
velocity plane with respect to the LSR 
(e.g. Oppenheimer et al. 2001; Reid et 
al. 2001). The squares (triangles) indi- 
cates the WDs with velocities in the 
Galactic plane (z=0 and 1^=0) larger 
(smaller) than 94 kms" 1 . The thick dot- 
ted and dashed circles correspond to the 
1 and 2— <r velocity dispersion contours 
for the disc and halo populations, re- 
spectively. The projection assumes that 
W=0. For the population with veloci- 
ties larger than 94 kms -1 , also the U 
and V velocities are shown under the as- 
sumption that v r = ± 30 kms - (dots 
with lines). The thin dotted ellipsoidal 
contours indicate the outermost orbital 
radii (-R max , in kpc) of WDs, whereas the 
paraboloidal contours indicate the inner- 
most orbital radii (-R m i n , in kpc). It is as- 
sumed that R=8 kpc at the Solar radius 
and v c =220 kms -1 . 



density. This should scale oc o"u an d it is found (Dehen 
Binney 1998b) that 



4" 



in Galactic Cartesian coordinates, where 



rt 



80 ±5 kms' 



(5) 



roughly consistent with a calculation using the disc density 
variation. Later in the paper (§5.1), we use this expression 
as a self-consistency check of our results. 

In addition, moving groups or non-axisymmetric den- 
sity variations (e.g. spiral structure) could be responsible 
for a second effect, a rotation of the velocity ellipsoid in the 
horizontal plane, known as vertex deviation. This effect is 
quite small for old stars (i.e. Ssl0°; Dehen & Binney 1998b) 
and we shall therefore ignore it. 

We notice here a key simplification that follows from 
the realization that we can observe stars moving with rel- 
ative azimuthal sky velocity V± out to a distance d max = 
V±//j,mi n . By contrast, the amplitude of the radial excursion 
is ^2QV/k 2 ~ 60d max . We can therefore ignore curvature of 
the stellar orbits and spatial gradients in the stellar density 
over the sample volume. A corollary is that we can also ig- 
nore a small rotational correction when relating the proper 
motion to the space velocity. 



nf is unit vector in the direction of the n-th star. The ve- 
locity that minimises the dispersion about its projection 
onto the measured tangential velocity v$ is easily shown 
to be [< A >~ 1 ]tj < v >j- Carrying out this operation 
for the full 99 WD sample gives a mean velocity, essen- 
tially the asymmetric drift velocity relative to the LSR, of 
(U, V, W) = (-15, -78, +9) ±12 km s -1 . This result is fairly 
robust in the sense that it is not seriously changed if we re- 
move a few very high velocities whose reality could be sus- 
pect. The sample has no significant motion either radially or 
perpendicular to the plane of the Galaxy, but clearly shows 
the asymmetric drift in V . Such an asymmetric drift, accord- 
ing to Eqn.fel), would naively imply a velocity dispersion for 
the complete sample of au ~ 79 kms -1 . The absence of a 
velocity drift in the U and W directions also supports the 
notion that this WD population is near equilibrium. 

The values of the velocity dispersion and asymmetric 
drift of a kinematically unbiased sample, however, will be 
lower than found above, because a proper-motion limited 
sample will preferentially bias toward the high- velocity WDs 
(see Sect. 5). 



4.3 Statistical Approach 

There is a standard procedure for finding the mean ve- 
locity of a sample of stars with known proper motions 
and distances, even when their distribution on the sky is 
far from isotropic (BM98). We define a projection matrix 



4.4 Galactic Orbits 

Because WDs can move distances from the sun that are at 
least ~2 orders of magnitude larger than their observed dis- 
tances (§4.2), the majority originate far from the LSR. To 
quantify this, we have calculated the inner and outermost 
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Figure 5. Examples of orbits for 
two WDs drawn from the sample, (a) 
High- velocity WD constrained to move 
in the Galactic plane with (U, V, W) = 
(-264,-169,0) kms" 1 . (b) Slow 
disk WD dwarf with (U,V,W) = 
(49, -21,0) kms" 1 . (c) The same WD 
as in (a) except that now we have as- 
signed the velocities on the assump- 
tion that the WD is drawn at ran- 
dom from the high- velocity component 
of the bimodal velocity distribution 
(§4.5). A projection onto the merid- 
ional plane is shown, (d) The same 
WD as in (b) with velocities drawn 
from the slow component. 



radii (r m i n and r max ) of the white-dwarf orbits in the simpli- 
fied Galactic potential, Eqn.(H), using the velocities U and 
V, and explicitly assuming W=0 (i.e. the WDs move only in 
the Galactic plane). The assumption H / =0 leads to a tight 
scatter of the WDs around the relation (U 2 + V 2 ) = (v 2 +v 2 ), 
whereas the assumption v r —0 leads to results are similar 
within a few percent. Statistically, these two different as- 
sumptions are therefore expected to lead to only minor dif- 
ferences. 

The results, shown in Fig. 4, indicate that as many as 
98, 72, 29 and 14 of the 99 WDs make excursions more than 
1, 3, 5 and 7 kpc, respectively. None of the WDs have both 
the inner and outermost radii inside the volume probed by 
the most distant observed WD (at ~150 pc). Hence, even un- 
der the restrictive assumption that these WDs do not move 
out of the Galactic plane, the majority are not of local ori- 
gin, but come from a volume at is at least 3 to 4 orders of 
magnitude larger than that probed by the survey. 

To illustrate this more graphically, we have plotted the 
orbits of two white dwarfs in the sample in Fig. 5, assuming 
firstly that the WDs lies in a thin, equatorial disk, i.e. W = 
0, and secondly that they are drawn from the bi-Maxwellian 
velocity distribution discussed below (§4.6). One WD has a 
high velocity and moves on a nearly radial orbit, the other 
moves on a nearly circular orbit. It is clear that the epicyclic 
approximation is unlikely to provide an accurate description 
for the orbits of most high velocity WDs. 



4.5 Disc Models 

In general the disc of our Galaxy is thought to be composed 
of two major dynamically distrinct components; the thin 
disc and a thick stellar disc, which supposedly comprises 
several percent of the stellar disc mass (e.g. Gilmore & Reid 
1983; Gilmore, Wyse, & Kuijken 1989). 

Whereas the thin disc has a relatively low velocity dis- 
persion of cri4/~20 km s -1 in the vertical direction and is only 
several hundred pc thick, the thick disc has velocity disper- 
sion that is <7w~45 kms - and has an asymmetric drift 
of i> a ~40 kms - (e.g. Gilmore et al. 1989; Chiba & Beers 
2000). If discs are treated as massive self-gravitating slabs, 
such that a^y ~ 2-kGTjZo, thick-disc stars typically move 
zo~l kpc out of the plane for the appropriate Galactic disc 
surface densities E, consistent with observations (e.g. Reid 
& Majewski 1993; Robins et al. 1996; Ojha 2001). 

Maybe surprisingly, recently analyzed observations sug- 
gest (i) that WDs which kinematically appear to belong to 
the thin disc, have a scale height that might be twice that of 
the traditional thin disc and (ii) that the WD sky-density to- 
wards to North Galactic Pole is an order of magnitude higher 
than found in previous surveys (Majewski & Siegel 2001). 
In comparison, about ~50% of the high proper-motion WDs 
move at least 4 kpc from the LSR (§4.4). Hence, WDs with 
velocities several times the velocity dispersion of the thick 
disc with respect to the LSR could, in principle, move many 
kpc out of the Galactic plane as well, if their velocity dis- 
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tribution function is not extremely flattened in the vertical 
direction. 

WDs observed in the Solar neighborhood are therefore 
a mixed population from the thin disc and the thick disc. 
Because of the proper-motion cut, we expect only a few WDs 
from the thin disc in the sample. If the halo contains a sig- 
nificant population of WDs, it will also contribute to the 
local white-dwarf density. Each of these populations have 
their own kinematic properties, although there could be a 
continuous transition, for example, from the thick disc to 
the halo (e.g. Gilmore et al. 1989). 

Whereas, epicycle theory (§4.2) indicates that excur- 
sions in the radial and tangential direction are related 
through conservations of angular momentum, unfortunately, 
we have little information about the velocities perpendicular 
to the plane. However, if most WDs are born in the disc, one 
expects that scattering processes in the plane of the galaxy 
result in comparable radial, azimuthal and vertical veloci- 
ties. This is observationally borne out by the fact that the 
velocity ellipsoid is isotropic within a factor ~2 (§4.2) and 
not extremely flattened. Similarly, if WDs occupy the halo 
as the result of a merger and/or violent relaxation process, 
during Galaxy formation, or if WDs are ejected into the halo 
from the disc (see §7.2) the velocity ellipsoid is also expected 
to be nearly isotropic. Theoretically nor observationally is 
there any evidence or indication that a velocity distribution 
can be extremely flattened in the vertical direction. 

If one accepts this conclusion, it is difficult to escape the 
consequence that the majority of the observed WDs with 
large U and V velocities also have large W velocities and 
make comparably large vertical excursions. Even, if we con- 
servatively assume that the vertical excursions are only half 
that in the radial and tangential direction, one finds that 
~50% of these 'local' WDs in the sample are expected to 
move ~ 2 kpc out of the plane, twice the thick disk scale 
height of zo~l kpc, although this conclusion could be modi- 
fied if a proper potential model for the stellar disc and bulge 
is included (e.g. §7.2). In the following subsection we will ex- 
amine this in more detail, in the context of a local velocity 
distribution function and show more generally how these re- 
sults can be translated into a global halo density of WDs. 



4.6 The Schwarzschild Velocity Distribution 

The small sample of WDs and lack of radial velocities does 
not allow us to derive the phase-space density of this pop- 
ulation, F(x, v), where x and v are position and velocity 
vectors, respectively (e.g. BT87). We must therefore make 
some simplifying assumptions: (i) the phase-space density 
is locally separable, such that F(x,v) = n(x) ■ f(v), (ii) 
the velocity distribution function (VDF) with respect to 
an observer at the LSR can be described by a normalized 
Schwarzschild distribution function (SVDF) 



fs(U,V,W) = 



exp(-u 2 /2) 
(27r) 3 / 2 au&vo'w 



(6) 



where v—(U/a v , (V + v a )/a v , W/a w ) and (iii) n(x) is con- 
stant in the volume probed by the survey. The adopted co- 
ordinates system is shown in Fig. 2. Both the velocity vector 
(U,V,W), with respect to an observer at the LSR, and the 
velocity dispersions au,(Tv and aw are defined along the 



principle axes of the Galactic coordinate system. In the re- 
mainder of this paper, we adopt a Bimodal Schwarzschild 
Velocity Distribution Function (BSVDF) model 



f B (U,V,W) = 



r exp(-tf? d /2) 



(1 + r)(27r) 3 / 2 o-[ /itd (Tv,tdO- W ',td 
exp(-v£/2) 



— (7) 

(l+r)(27r) 3 / 2 cra, h (Tv,i 1 aw,h' V ' 

where the definitions are as in Eqn.dq) and r — n t( j/nh is the 
number-density ratio of WDs in the thick disc and halo. The 
velocity dispersions, asymmetric drifts and number-density 
ratio are free parameters. We emphasize that the adopted 
BSVDF model is not a unique choice and the real VDF 
could be a continuous superposition of S VDFs or have a form 
different from Maxwellian. Experiments with more general 
forms of the VDF (power laws), however, give essentially 
similar results. In addition, the SVDF is locally a solution 
of the collisionless Boltzmann equation and also internally 
consistent with epicyclic theory (BT87). The same is true 
for a superposition of SVDFs. 



5 MAXIMUM LIKELIHOOD ANALYSIS 

To calculate the probability of a set of model parameters, 
M={au tl j l ,o'u,h, v a,d,Va,h,r} , given the set of constraints Di = 
{fj,i,pa,i}, i.e. the proper motion /i and position angle pa from 
the i-th white dwarf, we use Bayes' theorem 



Pi(M\Di 



P(M)P(Dj\M) 
P{Di) 



(8) 



We assume that P(M) is a constant and set it equal to unity, 
although this might not the case. In particular, the efficiency 
of detecting WDs might be a function of its proper motion 
and position on the sky (Oppenheimer, private communica- 
tion). If the values of P(Di\M) are statistically independent, 
the log-likelihood sample becomes 



C(M\D)= J2 lnP(A|Af), 



(9) 



where D = {Di} and A\vd is the number of WDs in the 
sample. The probability P(Di\M) for each WD in a fully 
proper-motion limited is given by the function 



P(Di\M) = 



XT W(vi,Vb)f(v = pi+ Vrh) dv r 
J v d 2 v J^ W(vi,Vb)f(v = Pi + v r fi) dv r 



(10) 



where pi=(ve,Vb)i is the velocity vector on the sky in Galac- 
tic coordinates (see Eqn.l) and fi is the unit vector in the 
radial direction away from the observer (Graff and Gould, 
private communications). The weight function W(vi,vt) = 
(vf + vf) 3 ' 2 and V indicates that area in velocity space for 
which (vf + vf) > (/iminf;) 2 1 where n is the distance to the 
WD and p m i n = 0.33 arcsec yr _1 is the proper-motion lower 
limit of the survey. The weight function is proportional to 
the maximum volume (V max ; see Sect. 6) in which the WD 
could have been found and therefore biases a proper-motion 
limited survey towards high- velocity WDs. In case /i m i n = 0, 
the likelihood function reduces to that of a kinematically un- 
biased survey (e.g. Dehnen 1998). 

In practice, only ~90% of the WDs in the sample are 
proper- motion limited, whereas only ~10% of the observed 
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Table 1. The most likely BSVDF model for the full (D g6 ). The 
errors indicate the 90% confidence interval. 



Parameter 



£>96 



<T(7,d (kms -1 ) 


62 +8 
DZ -10 


cr^h (kms -1 ) 


180±S 


v a , d (kms -1 ) 


5o±i; 


v a . h (kms -1 ) 


1761™ 


r = "disc /"halo 


16.0±»;° 



WDs would first drop out of the sample because of the sur- 
vey's magnitude limit (Fig. 7; see Sect. 6 for a more detailed 
discussion). These WDs are also those with highest observed 
velocities (Fig. 7). If we assume that these WDs are also 
proper-motion limited, we would overestimate the volume 
(i.e. W(vi,Vb)) in which these WDs could have been de- 
tected and thus bias the results towards lower velocities and 
lower halo densities (Gould, private communications), the 
latter by about 40-50%. Of the ten WDs that are magni- 
tude limited, eight have nearly identical values of V ma x (and 
therefore absolute magnitude), whereas the other two have 
much smaller values. Conservatively, we assume that all ten 
WDs have similar and constant values of Vinax, which might 
slightly underestimate the halo density of WDs. In prac- 
tice this effect is only ~10% and we therefore neglect it in 
light of all the uncertainties. It is important to remember, 
however, that strictly speaking Eq.(hOJ) is only valid for com- 
pletely proper-motion limited surveys. The notion that this 
is the case for ~90% of the WDs allows us to make the 
simple adjustment above without explicit knowledge of the 
white-dwarf luminosity function, which is roughly speaking 
assumed to be a delta function. 

For each WD, i, we now know the velocity vector p; 
(see Fig.3). We then evaluate P(Di\M) using Eqn.([uJ and 
the transformation between (vi,Vb, V r ) and (U, V, W), as dis- 
cussed in Sect. 3 and from those we can evaluate C(M\D), 
given a set of model parameters M. 



5.1 Results from the BSVDF model 

We maximize the likelihood C(M\D), given the BSVDF 
model introduced in §4.6 and the proper-motion data set 
from the 99 observed WDs. Before the optimisation, we first 
remove three of the WDs that contribute most (~80%) to 
the WD density determined from the complete sample. A 
justification for this is given in Sect. 6. The removal changes 
the results only a few percent compared to the dataset 
of 99 WDs. The resulting dataset from 96 WDs is des- 
ignated I?96. We optimize by varying all five parameters 
M={tJu,d,o'u,h,Va,d,Va,h,r}. For definiteness we assume that 
the velocity ellipsoid has ratios au'-av'-aw~^'-'2/3:l/2, close 
to what is observed for the smooth background of stars in 
the local neighborhood (e.g. Dehnen 1998; see §4.2). The 
resulting model parameters are listed in Table 1. The er- 
rors indicate the 90% confidence level, determined by re- 
optimizing, whereby we allow all parameters to vary, but 
keep the parameter of interest fixed. 

Fig. 6 shows the likelihood contours of thick-disc and 
halo velocity dispersions around the most likely models in 
Table 1. We maximize the likelihood for each set of velocity 



o 
o 



o 
o 

CO 



R o 
M o 

w C\2 
o 

b 

o 

o 




40 



60 



80 



<7 



disc 



(km s ) 



Figure 6. Likelihood contour plot of tJdisc versus a"halo for the 
full data set (96 WDs). The dashed line indicates a"di a c = "halo- 



dispersions. In the case Dge, one can conclude that VDF is 
bimodal at the more than 99% C.L. 

Due to the position of the WDs close to the South 
Galactic Pole, only poor constraints are obtained on the 
aw- Our particular choice of the ratio au '■ ffw—1 ■ 1/2 
could potentially underestimate the density of WDs in the 
halo. By varying the ratio (aw /&u) and re-optimizing the 
likelihood, we find (aw /au) < 1.1 at the 90% C.L., with no 
sensible lower limit. This exemplifies the poor constraints on 
velocities perpendicular to the Galactic plane. Within this 
range the parameters in Table 1 vary by only a few percent. 

Note also that the maximum-likelihood estimate of the 
thick-disc asymmetric drift, w a ,d=50 kms -1 , agrees well 
with that estimated from the asymmetric-drift relation in 
Eqn.(5), which gives 48 kms -1 . If we make the velocity vec- 
tor isotropic, for example, this is no longer the case and 
the resulting asymmetric drift and radial velocity disper- 
sion of the thick-disc population (i.e. crry=51 kms -1 and 
u a =49 kms -1 ) are inconsistent with the empirical relation 
in Eqn.(5). The likelihood of the model also decreases by 
a six orders of magnitude. This might be regarded as a 
consistency check that the sample of thick-disc WDs and 
the likelihood results are in agreement with the kinematics 
of other local stellar populations, whereas this is not the 
case for an isotropic velocity distribution. The flattening of 
the velocity ellipsoid in the V direction also naturally ex- 
plains the apparent scarcity of WDs in the regions around 
(U,V) = (0,-220) kms -1 and V < -220 kms -1 , noticed 
by Reid et al. (2001). 
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Figure 7. The ratio Vmax/Knax 5 as function of observed ve- 
locity on the sky (iijSVg + v^). The ten WDs (squares) with 
Knax/Vmax g > 1 are magnitude-limited in the sense that they 
would still be in the proper motion sample if they were at a 
distance such that their apparent magnitude equaled the survey 
limit. The remaining 89 WDs (circles) are proper-motion lim- 
ited. Note that most of the magnitude-limited WDs have high 
velocities and low luminosities. A linear fit to sample shows a de- 
pendence XV >■? , as expected for a proper-motion limited sample 
with no correlation between white-dwarf luminosity and velocity. 



Figure 8. The correlation between velocity on the sky (v±) and 
magnitude R59F. The lower dotted line results from the rela- 
tions V^ x = KjJax, /i m i n =0.33 arcsec yr _1 and a constant abso- 
lute WD magnitude < Mn^gp >=14.3 mag, i.e. the average of 
the sample of WDs. Below this line, WDs are either too faint or 
slow to be detected. Some WDs are brighter than Mr$qp=14.3 
mag and lie below the line, but in general most stars follow this 
lower limit. The upper dotted line indicates the same, but for 
A'min=l arcsec yr _1 . The local WD density is dominated by the 
3 WDs (stars) that lie isolated from the bulk of the population, 
and which presumably belong to the thin-disc (see text). 



6 THE LOCAL WHITE DWARF DENSITIES 

We must now turn to estimating the local space density 
which is needed to normalize the stellar distribution func- 
tion. Our approach follows that of Oppenheimer et al (2001) 
by taking into account the proper motion selection, as well as 
the flux selection. We continue to ignore errors, which can 
introduce significant biases, and suppose that a star with 
proper motion [i, magnitude R59F and distance d, could 
have been seen out to a distance d max = min[dmax,dmax] 
where the proper-motion distance limit is d^ ax = Vj_//x m in, 
with /i m i n =0.33 arcsec yr~ , and the magnitude distance 
limit is given by log((Cax) = 1 + (19.8 - R59F)/5 and still 
been included in the sample. 

The maximum volumes in which the WD could have 



been found are then V, 



]ijr, 



(a/3) x (d^/ax) 3 , where O a is 
solid angle covered by the survey. In Fig.7, we plot the ra- 
tio VJ^x/V^ax versus observed velocity (v±). One observes 
that 89 out of the sample of 99 WDs are proper motion- 
limited; there is a conspicuous absence of faint WD in the 
sample relative to a true magnitude-limited sample. Note 
that the absence of faint WDs is even more prominent for the 
slowly moving stars (Fig. 8). There are two possible explana- 
tions. Either the photometry is inaccurate and the selection 
is quite incomplete close to the survey limit R59F = 19.8. 
Alternatively, the density of old, low luminosity WDs is so 
small that there are no more than a few that are close enough 
for us to see them move. For this second explanation to be 



viable, the logarithmic slope of the luminosity function $(£) 
must be greater (and in practice much greater) than —5/2. 
In order to distinguish between these two explanations we 
have evaluated < V/Vmax >~< (d/dmax) 3 > for the whole 
sample. We find a value < V/Vm&x >— 0.41. This is quite 
close to the expect value of 0.5 that we suspect that although 
incompleteness is present, the distribution in d — v — L pa- 
rameter space reflects the whole population, where L is the 
luminosity of the WD. 

If we ignore possible photometric incompleteness, our 
estimator of the total space density is 



total 
"0,WD 



3/nVs 

47T 



E.-3 
max,i 

i=l 



(ii) 



where /n = 0.12 is the fraction of the sky surveyed and 
fs — 0.8 is the estimated spectroscopic incompleteness of 
the 99 WD sample. The result is n °wrj = 0.0096 pc~ 3 . The 
result is quite robust with respect to changes in the limit- 
ing magnitude although 30% of the density is contributed 
by a single WD, and 80% by three WDs (i.e those with 
R59F < 14 mag; Fig. 8). This estimator refers to stars whose 
luminosities and transverse velocities are actually included 
in the sample although the variance can be quite large if 
there are significant contributions from WD near the bound- 
aries of this distribution. From Fig.l it is apparent that this 
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range translates into absolute magnitudes, 16 ^ Mr^f ^19 
and ages in the range of about 1 — 10 Gyr. 

The three brightest WDs that dominate the density es- 
timate are most likely part of the thin-disc population. Be- 
cause this component was not included in the VDF model 
used in the likelihood analysis, we have to remove these ob- 
jects from the density estimate to avoid a severe overesti- 
mate of rig ^ D , the normalisation of the local phase-space 
density, F B (U, V, W) = n* d +^ x f B (U, V, W), of the thick- 
disc and halo WDs. We therefore remove all WDs with ob- 
served velocities v± < 30 kms -1 from the sample (three 
in total), which should remove all thin-disc WDs. Their re- 
moval from the sample has a negligible effect on the like- 
lihood analysis (see §5.1). We subsequently find nf^ 1 ^ = 
(1.8 ± 0.3) x 10~ 3 pc~ 3 (1-c), which then represents the 
local density of thick-disk and halo white dwarfs that have 
v± > 30 kms . 

To normalise Fb{U,V,W) we calculate the fraction, 
/>30, of the VDF /b(U, V,W) for which the observed ve- 
locity v± > 30 kms -1 . The correct normalisation of the 
local density of thick-disk and halo WDs then becomes 



td+h 



-td+h 

l 0,WD 



//> 



Given the BSVDF model listed in Ta- 
ble 1, for the full data-set (Dge), we find < />3o >= 0.946, 
hence 



™o d WD = (1-9 ± 0-5) x 10~ 3 pc~ 3 



(90% C.L. 



Finally, the local halo density of WDs becomes 

l.llo'.T X 10 " PC 



h 

%,WD 



(1 + r) 



- ' „ in-4 „„ 3 ^ 9Q% CL 



with r=16.0+io5 (Table 1). This density estimate for the 
halo white dwarfs is about half that found by Oppenheimer 
et al. (2001). For completeness, we quote a thick-disc white- 
dwarf density of 







Figure 9. Three dimensional WD velocity distributions at two 
locations. The top row corresponds to radius r = 0.7r c = 5.6 kpc 
and colatitudc 9 = 60° relative to the Galactic center and the 
bottom row to r = 1.4 r c = 11.2 kpc, 9 = 60°. The projections on 
the left are of f(v r , 0, v^) and those on the right are of f(v r , vg, 0). 
The velocity scale v c = 220 km s — 1 is marked. The four contour 
levels are at 0.01, 0.03, 0.1, 0.3 times the local peak halo WD dis- 
tribution function, /o = 1.5 X 10 -11 pc -3 (kms -1 ) -3 . The shaded 
areas are excluded volumes in velocity space where, under the 
dynamical assumptions discussed in the text, local measurements 
of the WD distribution cannot inform us about the halo distribu- 
tion. 



td 
n 0,WD — 



(1 + r) 



= (1.8 ± 0.5) x 10~ 3 pc~ 3 (90% C.L.), 



which corresponds within 1.5 <r to the thick-disc number den- 
sity rio d WD — (3-3 ± 1.2) x 10~ 3 pc~ 3 that was estimated by 
Reid et al. (2001) on the basis of only seven WDs within a 
sphere of 8 pc from the Sun. Note that our estimate takes 
the contribution from WDs outside the 94 kms -1 cut in 
to account. We therefore self-consistently recover the local 
thick-disc density within acceptable errors, from the sam- 
ple of high-proper motion WDs discovered by Oppenheimer 
et al. (2001). We therefore do not expect this sample to be 
significantly incomplete or biased in a way unknown to us. 
Concluding, the likelihood result that there are two 
kinematically distinct populations, at the 99% C.L, gives 
us confidence that there is indeed a contribution of ~6% of 
the thick-disk WD density or ~0.8% of the local halo den- 
sity, by a population of WDs that exhibits kinematic prop- 
erties that are not dissimilar to that expected from a nearly 
pressure-supported halo population. 



6.1 Local Stellar and Halo Densities 

To place the results obtained above in perspective, we now 
turn to a comparison of the local density of WDs that wc 
derived with other local densities. 

Summarising, the local mass densities of WDs, derived 



in this paper, were Po.wd 
6 x 10~ 3 M Q pc~ 3 



0.7 x 10" 



1.1 x 10" 3 and 



respectively, for the halo, thick-disc 
and total mass density (including thin disc), if we convert 
number densities to densities using an average WD mass of 
0.6 Mq. The local WD density derived from direct obser- 
vations is (3.7 ± 0.5) x 10~ 3 M Q pc~ 3 (Holberg, Oswalt & 
Sion 2001). These direct observations compare well with our 
total WD density estimate. Because of the very small vol- 
ume probed by the proper-motion limited survey of Oppen- 
heimer et al. (2001) for the lowest velocity thin-disc WDs, 
the error on our estimate is very large. We found ~3 plausi- 
ble thin-disc WDs (Sect. 6; Fig. 8). If we had found one WD 
less, with comparable or lower velocities, the total density 
estimate could easily have halted. The low number statistics 
of thin-disc WDs, which dominate the local WD density, 
makes the integrated density estimate of local WDs very 
uncertain. This is much less of a problem for both the thick- 
disc and halo density estimates, because in both cases the 
probed volume is significantly larger (2-3 orders of magni- 
tude) and consequently also the number of WDs from both 
populations. 

The local stellar density is around 0.06 — 0.1 Mq pc~ 3 
(e.g. Pham 1997; Creze et al. 1998; Holmberg & Flynn 2000). 
The mass density ratio of WDs to stars is therefore about 
4-6%. This fraction is 2-3 times less than the ~13% es- 
timated by Hansen (2001). We find ~15% from a more so- 
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phisticated population synthesis model that self-consistently 
takes metallicity-dependent yields, gas in- and outflow, stel- 
lar lifetimes and delayed mixing into account, under the 
same assumption that the stellar IMF is Salpeter from 0.1 
to 100 Mq . The result is almost independent of any assump- 
tion, except for the shape and lower mass limit of the IMF. 
We return to this apparent discrepancy in §7.2. 

The local halo density in dark matter is estimated to be 
(8±3) x 10" 3 M pc" 3 (Gates et al. 1995; where we take the 
result that does not include the microlensing constraints). 
This result agrees with those determined by other authors. 
Our local halo WD density estimate would therefore cor- 
respond to ~0.8% of the local halo mass density, slightly 
higher than that estimated by Oppenheimer et al. (2001), 
for the reasons outlined previously. 

The local density in hydrogen-burning halo stars is es- 
timated to be (4 ± 1) x 10~ 5 M Q pc~ 3 , whereas the WD 
density estimated from this population, based on popula- 
tion synthesis models, is 1.3 x 10~ 5 M pc -3 (Gould et al. 
1998). Our density estimate of WDs is therefore ~5 times 
larger than that previously estimated and is comparable in 
mass to the local stellar halo density. Such a high density in 
halo WDs demands a non-standard explanation, other than 
that based on simple population synthesis models and/or a 
standard IMF. We discuss this in more detail in §7.2. 

First, however, we examine the contribution of these 
high- velocity WDs to the non-local halo mass density, which 
allows us to estimate their total mass contribution to the 
Galaxy, under some simplifying assumptions. 



7 THE HALO WHITE DWARF DENSITY 

Combining the results of Sections 5 and 6 we now have a 
normalised halo distribution function in the form 



F h (x ,v) = n , WD x 



exp(-« 2 /2) 
( 2-7T ) 3 / 2 au ay aw 



(12) 



evaluated at the solar radius (r c — 8 kpc). What does this 
tell us about the total white-dwarf density throughout the 
halo? The connection is provided by the collisionless Boltz- 
mann (or Vlasov) equation which states that the distribu- 
tion function does not vary along a dynamical trajectory. 
Therefore, if we knew the halo potential <I>(r) as well as the 
local velocity distribution function, we would know the halo 
distribution function at radius r and colatitude 9 for all ve- 
locities that connect to the solar radius. (We assume that 
the distribution functions is symmetric.) 

We now make a further simplification and assume that 
the Galactic potential, &(f)=&(\r\), is spherically symmet- 
ric. Specifically, we adopt Eqn.(H) although our results are 
not particularly sensitive to this choice. In making this as- 
sumption, we are ignoring the effects of the disk matter 
which will retard the WDs as they move to high latitude, 
an effect that is relatively small for the high velocity, halo 
stars of interest. However, this is a bad approximation for 
the two disc components. 

A WD on a high latitude orbit will oscillate in radius be- 
tween a minimum and a maximum radius determined by en- 
ergy and angular momentum conservation. As the orbit does 
not close, the ap- and perigalacticons (r maXjm i n ) will eventu- 
ally become arbitrarily close to the Galactic disk. Therefore 



as we are assuming axisymmetry, (and the orbital preces- 
sion that would be present in a more realistic non-spherical 
potential will obviate the need for this assumption), we are 
locally sampling the halo velocity distribution at (r, 9) for 
halo velocities (v r ,Vg,V<f,) for which r max < r c < r m i n . For 
exterior radii, the region of velocity space that is not sam- 
pled is a hyperboloid of revolution of one sheet, 



(r 2 /r 2 - l)(vg + vl) - V 2 > 2v 2 ln(r/r c ). 



(13) 



For interior radii, the excluded volume is an oblate spheroid 



+ (1 - r 2 /rl)(v 2 e + vl) < 2v 2 ln(r c /r). 



(14) 



Two examples of exluded regions are shown in Fig. 9. 

The connection between the solar neighborhood and the 
halo is effected using the constants of the motion (e.g. May 
& Binney 1986), for which a convenient choice is the energy, 
squared angular momentum and the angular momentum re- 
solved along the symmetry axis: 



E 



^(U 2 + (v c + V) 2 + W 2 ) + $(r c ) 



= r 2 c [(v c + V) 2 +W 2 ] 
= r c (v c + V). 



(15) 



This enables us to compute F(E, L 2 ,L Z ) in the halo at (r, 
within the sampled volume of velocity space, using 



v 2 - 2E - 2$(r) - L 2 /r 2 
L 2 /r 2 -L 2 Jr 2 sin 2 8 
L z /r sin 9. 



2 

I'D 



(16) 



Making these substitutions, we discover that the velocity 
distribution function in the halo at (r, 8) is also a Gaussian, 



F h (r, 8;v r ,v^,Vg) 



"o.wd 



(2ir) 3 / 2 auavaw 



x exp 



A$ 



x exp 



2a% 



x exp 



-vl 
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x exp 



_£ 
2 

x exp 



\-x z 



+ ■ 



i 

+ - 



u u u v 

V${v c — Va)xsm.t 



(y c - v a y 

2<4 



(17) 



where A"I> = $(r) — &(r c ) is the difference in the po- 
tential between the point (r, 9) in the halo and the Solar 
neighborhood. Our particular choice for the potential gives 
A$ = v 2 ln(a;), with x — r/r c . Although this choice for the 
potential is particularly convenient, the above VDF for the 
halo can be used with any spherical symmetric potential 
model and locally constrained Schwarzschild velocity distri- 
bution function. Now in order to compute the density, we 
must make an assumption about the velocity distribution 
function within the excluded volume. The natural assump- 
tion is that it is the same Gaussian function as within the 
connected volume. We emphasize that this is probably a 
pretty good assumption when the excluded region of ve- 
locity space is small, though quite suspect when it is not. 
Performing the elementary, if tedious, integrals, we obtain 



n(r,9) 



n h (r c ,n/2)^_ (Va/au) 2 



ayaw 



-,-1/2 



+ 
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Figure 10. Left: The halo WD density derived from the most likely BSVDF model (e.g. Table 1) and Eqn.(hq). Contours indicate 
factors of two increase (solid lines) or decrease (dashed lines) in density compared to pi} WD = 0.7 X 10 — 4 M© pc ~at the LSR at radius 
-R=8 kpc. Right: The same for the thick-disc WD density with PQ d WD = 1.1 X 10 -3 Mg pc -3 . Note that the true thick-disc WD density 
distribution should be more affected by the non-spherical disc potential than the halo population, and the computed thick-disc density 
should only be regarded as indicative. 



+ 



+ ■ 



-1/2 



(18) 



exp 



(V a - VcY 



(1-^) , x 2 cos 2 



2ol 



(1 — x 2 ) , x 2 sin 2 
1 T 1 — 

"u a v 



where the factor (r/r c )~^ v "' au ' comes from the poten- 
tial, i.e. exp[— A$/cry]=exp[— v^/afj ln(r/r c )]. This proce- 
dure can only be self consistent for au < Max[ov, aw] when 
r < {1 — (Jy/Max[cry, crj^]}^ 1 ' 2 , otherwise the computed 
density is unbounded. (Indeed there appears to be a physical 
limit on the vertical velocity dispersion of a halo population 
within the Galactic plane.) 

The density distributions of the halo and thick disc of 
the most likely model are plotted in Fig. 10. The halo den- 
sity distribution can be seen to be roughly in the form of 
an oblate flattened spheroidal distribution with axis ratio, 
ijwd ~ 0.9 at the solar radius. We expect q is decrease 
slightly if a more realistic flattened potential is used, in- 
cluding constributions from the bulge and disc (§7.2). If we 
assume that the average WD mass is 0.6 Mg, then the mass 
contained within a shell 0.3-3 r c is ~1.5xl0 9 Mq. The to- 
tal WD mass inside 50 kpc is ~2.6xl0 9 M© for our most 
likely model. The radial variation of density is n a r - 30 
as argued on the basis of the fluid model and the mass per 
octave varies only slowly. This result shows that the mass 
fraction of WDs decreases as oc r _1 '° , in the particular po- 
tential that we adopted. Of course the bulge and the disk 
must certainly be included within ~ 3 kpc and we do not 



expect the circular velocity to be remain roughly constant 
beyond ~ 30 kpc and our model is quite primitive, relative 
to what could be created with a more extensive data set and 
a more sophisticated form for $(r). 

Even so, if we were to believe the simple potential 
model, the total WD mass corresponds to ~0.4% of the to- 
tal mass enclosed by 50 kpc, Hence, we find S7wd ~ 10~ 4 , if 
fihaio(< 50 kpc) ~ 0.03 (e.g. Bahcall et al. 2000). This corre- 
sponds to ~0.4% of the total baryon budget in the universe. 
Clearly the halo WDs are dynamically unimportant. 

7.1 The Luminosity of Halo White Dwarfs 

A halo mass of ~2.6xl0 9 Mq in WDs corresponds to 
-/Vwd ~4.3xl0 9 WDs in total. The expected luminosity of 
this population in iZ-band is 

Rwu « 5 log(dp C ) + (Mr) ~ (5/2) log iV WD - 5 (mag), 

where (Mr) is the average absolute _R-band magnitude of 
the sample of WDs. If we take (Mr) ~ 14.3 from the sample 
of observed WDs by Oppenheimer et al. (2001), this reduces 
to -Rwd ~ 51og(d pc ) — 14.8. Similarly, if the luminosity of 
halo stars is Lr ~ 4x 10 7 L© (BM98; conveniently assuming 
Solar colors for the halo stars) and Mh,q=4.3, one finds 
Rst&is ~ 51og(d pc ) — 20. If these values are typical for L* 
galaxies, then the integrated luminosity of halo WDs in R- 
band will be ~5 mag fainter that that of halo stars, both of 
which are of course significantly fainter than the stellar disc. 
We conclude that it will be exceedingly difficult to observed 
such a faint halo around external galaxies. This is made even 
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harder, because this emission is spread over a much larger 
solid angle on the sky than for example the stellar disc. 

7.2 The Origin of Halo White Dwarfs 

We now turn the crucial question: if the density of halo WDs 
is as high as the observations, taken at face value, and our 
subsequent calculations indicate, where did all these WDs 
come from? The apparent similarity in magnitude, color and 
age between the high-velocity halo and low-velocity thick- 
disc WDs (Fig.l; Hansen 2001), suggests that the halo WDs 
originate in the disc and are subsequently ejected into the 
halo, through a process yet unknown to us. Their disc-origin 
is supported by the low density of halo stars (e.g. Gould et 
al. 1998). If halo stars and WDs come from the same (halo) 
population, there should be significantly more stars than 
WDs in the halo for any reasonable IMF. 

One possibility is that binary or triple star interactions 
eject stars into the halo with high velocities (socalled run- 
away stars), although not preferentially WDs. Alternatively, 
it is possible that a large population of WDs were born in 
globular clusters on radial orbits or satellite galaxies that 
have been destroyed by tidal forces over the past 10 Gyr, 
through which most of their stars end up in the bulge. Both 
scenarios are speculative and also do not satisfactory explain 
why the ratio of WDs to stars is so high in the halo. 

We propose another possible mechanism, in which most 
of the high-velocity WDs orginate from unstable triple 
stellar systems. Observations indicate that a high fraction 
(~30%) of binary stars in the Solar neighborhood have a 
companion and are in fact triples (e.g. Petri & Batten 1965; 
Batten 1967). If for example the third outer-orbit star, or 
one of the binary stars, evolves into a WD or goes supernova, 
or if matter can flow between the different stars (e.g. through 
Roche-lobe overflows and/or stellar winds) an initially sta- 
ble system can change its ratio of orbital periods and at 
some point become unstable, if the ratio of orbital periods 
(Pi/P s , the ratio of the largest over smallest orbital period) 
exceeds a critical value (e.g. Anosova & Orlov 1989; Kisel- 
eva, Eggleton, & Anosova 1994; Anosova, Colin, & Kiseleva 
1996) that depends on the mass ratio of the outer star and 
the inner binary. The outer-orbit star, often an evolved star 
itself, can then be ejected from the triple system with a high 
spatial velocity (e.g Worrall 1967), giving the inner compact 
binary system a "kick" in opposite direction, with a typical 
velocity observed for compact binaries in the disc (e.g. Iben 
& Tutukov 1999) and comparable to the orbital velocity of 
the compact binary. A second possibility is that the most 
massive star of the inner binary evolves into a WD, with- 
out the binary merging (Iben & Tutukov 1999), changing 
the orbital period and making the system unstable. Also in 
this case the WD can be ejected, leaving behind a recoiled 
binary. 

Besides explaining why even young WDs (few Gyr) are 
observed with high velocities, a consequence of such a mech- 
anism is that the total mass in high-velocity WDs (e.g. those 
in the halo), must be related to the integrated stellar mass 
in our Galaxy and the fraction of the stellar IMF in unsta- 
ble multiple (e.g. triple) systems. The estimated stellar mass 
in our Galaxy is roughly 6 x 10 10 M (BT87). If ~15% of 
the stellar mass equals the total galactic WD mass, based 
on simple population synthesis models (§6.1), the inferred 



Galactic WD mass should be ~ 9 x 10 9 M . In §6.1, we es- 
timated from observations that the local WD mass fraction 
was only ~5%, which then corresponds to ~ 3 x 10 9 Mq 
if this fraction holds throughout the disc. We estimated 
~ 2.6 x 10 9 M Q for the mass in halo WDs (Sect. 7), which 
then adds to a total WD mass of ~ 6 x 10 9 M . 

There are, however, a number of difficulties associated 
with the process described above. First, it implies that one 
out of every two WDs must be ejected from the disc, with 
an average rate of one WD every ~4 years, for an assumed 
Galactic age of ~ 10 Gyr. Second, it is not known what the 
birth rate and orbital parameters of triple (or higher-order) 
systems are. In particular, a large fraction of stellar systems 
could be intrisically unstable when they form. Such systems 
have been observed (e.g. Ressler & Barsony 2001), although 
their ultimate fate is unknown. The instability of triple sys- 
tems is particularly evident when comparing their orbital 
periods (e.g. Tokovinin 1997; Iben & Tutkov 1999), which 
shows a complete absence of triple systems with P\/P B £•> 2.5 
(i.e. the lowest theoretical limit for stability; e.g. Kiseleva 
1994), whereas there is no a priori reason to suspect that sys- 
tems with these orbital period ratios are not formed. A more 
typical ratio P\/P s where systems becomes unstable is ~4. 
Examining the distribution of Pi/P s in more detail reveals 
some clustering of triple stellar systems close to the stability 
limit for larger orbital periods (Tokovinin 1997; Iben & Tu- 
tukov 1999), although it is not clear whether that could be a 
selection effect. For the lower period systems this clustering 
is not the case. Third, because evolving triple stars have not 
been studied in as much detail as binaries, for example, the 
velocity of the ejected WDs is not well known in these pro- 
cesses, although it can be estimated from analytical models 
and numerical simulations. 

As a rough rule of thumb, stable triple systems can 
be separated into a close binary (masses mi and m.2) and a 
more distant companion (mass mz) with period ratio exceed- 
ing ~ 4 (see above). Let us consider two options. Firstly, sup- 
pose that mi is the most massive star and that it evolves to 
fill its Roche lobe and transfer mass, conservatively, onto iri2- 
This will initially shrink the orbit and increase the period 
ratio until m\ = rri2 when the separation is a m i n . Further 
mass transfer will then lead to an increase of the close binary 
orbit and a reduction of the period ratio oc m^m^ 3 '. As the 
close binary orbit expands, the proto-WD speed Vi will ini- 
tially increase oc m\m% and then decrease. For example, with 
mi = 2 Mq, m,2 = 1 M initially the close binary period 
could decrease by a factor ~ 3 by the time mi ~ 0.6 M 
and, provided that the companion's orbital radius jj 6o m m, 
the orbit can become unstable at this point. If mi is ejected 
by a slingshot effect involving ms, then a characteristic ve- 
locity at infinity is Vi ~ 200- (a m i n /4R ) -1 ' 2 . Higher speeds 
are possible with favorable orbits, especially retrograde com- 
panion orbits. Secondly, suppose that m-j, > mi + m,2. The 
mass transfer is less likely to be conservative although there 
is likely to be a net loss of orbital angular momentum asso- 
ciated with the loss of mass. The net result is that the outer 
binary orbit will initially shrink, again promoting instability. 

General analytical and numerical analyses of triple-star 
systems, randomly placed in a spherical volume, show that 
most (<;80%) systems are unstable and typically eject the 
lowest mass star with a median velocity comparable to their 
orbital velocity (e.g. Standish 1972; Monaghan 1976a,b; 
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Sterzik & Durisen 1998), as also indicated above. A more 
detailed analysis by Sterzik & Durisen (1998), where the 
ejection velocity distribution is weighted by the stellar mass 
function, shows that ~10 percent of the ejected stars (e.g. 
WDs) have velocities *J3 times the typical orbital velocity. 
Given their stellar mass function, the typical mass of an 
evolved triple-star systems is then ~5 Mq . As we saw previ- 
ously, if the orbital radius is of the order of the stellar radius 
(i.e. orbital periods of several days), escape velocities of sev- 
eral hundred kms -1 can be reached. For orbital periods of 
several hundred days (e.g. Tokovinin 1997), the typical or- 
bital velocity will still be ~80 kms -1 . WDs have to gain a 
velocity of ~ 100-150 kms -1 (in addition to their orbital ve- 
locity and velocity with respect to their LSR) to atain veloc- 
ities comparable to the observed high- velocity WDs. Hence, 
the proposed mechanism quite easily reaches the required 
velocities. Even in those systems with longer orbital periods 
(~1000 days), several tens of percent of the escaping stars 
will atain very high velocities (e.g. Sterzik & Durisen 1998). 

Additional support for the notion that these processes 
might be quite common, could be the good agreement which 
Valtonen (1998) recently found between the observed mass 
ratio of binaries and that determined from a statistical 
theory of three-body disruption (Monaghan 1976a,b). This 
might indicate that many stars are born in triple or higher- 
order systems, which become unstable and eject stars until 
only a stable binary is left. 

In order to explore this hypothesis further, we have con- 
structed a relatively simple model of the Galaxy as follows; 
we suppose that the mass comprises a spherical halo 



Ph(r) = 



(190 kms" 1 ) 2 3rl + r 2 



AtyG 



(r 2 +r 2)2 



(19) 



= 8 kpc, which generates a logarithmic potential, 
an 



with rh 
and a disc 

Pd(i?) = 60-e-- R/ ' ld <5(2) M Q pc~ 3 , 



(20) 



with fed = 3 kpc and a bulge, which we for simplicity flatten 
onto the Galactic plane, 

-3 



p b (R) = 190-e- R/hb 8{z) M 



Q pc 



(21) 



with /ib = 0.8 kpc. The combination of these potentials gives 
a nearly flat rotation curve beyond several kpc (e.g. Dehnen 
& Binney 1998a), with circular velocity v c w 220 kms~ . As 
the potential is no longer spherically symmetric, the total 
angular momentum is not an integral of motion. 

Given this model of the Galaxy, it is then a straigh- 
forward matter to follow the trajectories of observed WDs 
back in. In order to measure the contemporary WD source 
velocity distribution (per unit disk mass, per unit time, per 
unit speed), which we denote by g(v), we must make some 
assumptions. To keep matters simple, we suppose that g(v) 
is independent of radius and decreases with lookback time, 
t oc exp[— t/r], where r = 7.5Gyr as indicated by the data 
exhibited in Fig. 1. As individual orbits quickly sample a 
large range of disk radii it is not easy to derive the radial 
launch rate from the local WD distribution function. The 
best approach utilizing the present data set is to fit the halo 
distribution function from Table. 1. 

We use again Jeans' theorem to relate the observed 
VDF F h (xo,vo) to that at the source. Adopt one point in 
local velocity space and relate to g(v) through 
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where the sum takes account of m disk crossings along the 
orbit. The weights Wij are given by 



Wn = 



and n 



E(r-j) exp[-tj/r] 



4nv*Vj Z 



(23) 



and tj are evaluated at the j-th disk crossing. 
We next discretize g{vj) using 50 kms -1 bins and develop 
a least squares solution to Eq.(22) 



9 i = 






(24) 



where we sample the distribution function at N random 
points. 

The result indicates that the mean speed of the WD at 
launch is very high, w=250-300 km s _1 , with a birth rate 
of ~10~ 4 (Mq kms -1 Gyr) -1 at these velocities, consistent 
with ~10% of the stellar mass in the Galactic disc having 
been transformed in WDs. The low rotational velocity of 
the observed WDs requires them to originate in the inner 
2-3 kpc of the Galactic disk, so they loose most of their 
rotational velocity when climbing out of the Galactic poten- 
tial to the Solar neighborhood. The latter result is not un- 
expected, because most stellar mass (~75% for a disk scale 
length of h ~3 kpc) lies inside the solar radius and in addi- 
tion AM/Ar oc rexp(— r/h) peaks at r = ft. Hence a typical 
WD is expected to originate in the disc at r — h ~ 3 kpc. 

These calculations are consistent with, though do not 
require, that the WDs are a product of stellar evolution in 
the disk. Three lines of evidence support this conclusion, (i) 
The age distribution is just about what would be expected 
from a population of main sequence stars with a reasonable 
mass function from ~0.8 Mq to ~6 Mq. The stars are born 
at a rate consistent with stellar evolution. The most popu- 
lous stars are of mass ~ 1 — 2 Mq and take ~ 5 — 10 Gyr to 
evolve. Even if the Galaxy is ~ 14 Gyr old, we do not expect 
a large proportion of stars with ages <; 10 Gyr. A more care- 
ful, Galactic synthesis model is consistent with these state- 
ments, (ii) The velocities are isotropic with no discernible 
asymmetric drift. This is just what would be expected from 
a disk population that launches WDs in random directions. 
The median speed observed is ~200-300 km s _1 . (iii) The 
radial variation of white-dwarf sources appears to fall with 
radius in a similar fashion to the the stellar distribution. 

Of course, the situation is far more complex than this. 
However, the similarities between the characteristic escape 
speed and the halo velocity dispersion, together with the 
shortage of viable alternatives makes this a scenario that is 
worth investigating further. But it should again be empha- 
sized that a much larger sample of halo WDs is needed to 
place the above simple estimates on a more firm basis. 



8 MICROLENSING OPTICAL DEPTHS 

The next question to ask is whether the population of WDs 
in the halo can explain the observed microlensing events 
towards the Large Magellanic Cloud, with optical depth 
-tlmc = 1.2 x 10~ 7 (Alcock et al. 2000). Using the local den- 
sity and the global density profile of the population of WDs 
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in the halo, derived in the previous sections, it is possible to 
integrate over the density profile along the line-of-sight to- 
ward to LMC {£ = 280.5°, b = -32.9° and D = 50 kpc) and 
obtain the optical depth for any arbitrary set of parameters 
for the local VDF. 

For the most likely model, listed in Table 1 (full data 
set), we compute a microlensing optical depth of the halo 
WDs, twd ~ 1-3 x 10~ 9 . The optical depth of the thick 
disc is i~wd ~ 4 x 10 -9 . Hence, even in the most optimistic 
case, the integrated thick-disc plus halo WDs have an op- 
tical depth that is still 1-2 orders of magnitude below that 
observed. This difference can not be made up by a, most 
likely small, incompleteness of high- velocity WDs. 

If we scale the results in Section 6 and 7 to the LMC it- 
self, assuming (i) a similar shape of the velocity ellipsoid, (ii) 
velocities that scale with a factor (w^ MC /u^ al ) = (72/220) 
(i£ MC = 72 kms -1 ; e.g. Alves & Nelson 2000) and distances 
that scale as vj?/r=constant, and (iii) the ratio of total WD 
mass to stellar mass in the LMC is similar to that inferred 



in our Galaxy (~9% from Mttl 



5.3 xlO 9 M Q ;e.g. Alves 



& Nelson 2000), we find a self-lensing microlensing optical 
depth for the center of the LMC of r ~ 2 x 10~ 8 , lower than 
the results from the MACHO collaboration (Alcock et al. 
2000). We find that the results are quite sensitive to the as- 
sumped parameters, especially the velocity dispersion. The 
WD density is oc r _ '" c//CTC/ - ) and gets unbounded if the ve- 
locity dispersion is lower somewhat. This illustrates that a 
"shroud" of white dwarfs around the LMC (e.g. Evans & 
Kerins 2000), in principle, could account for the observed 
optical depth, although it does require that 50% of the halo 
WD mass is concentrated in the inner 0.7 kpc and rather 
fine-tuned model parameters. The latter could be overcome 
if the mass profile of the LMC has a finite core. 

A similarly rough scaling for the edge-on spiral gravi- 
tational lens B 1600+434 - where we assume 4 times more 
mass in WDs than in B1600+434 - suggests t W d ~ 10~ 3 
along the line-of-sight towards lensed image A, at 6 kpc 
above the galaxy plane. This optical depth appears insuf- 
ficient to explain the apparent radio-microlensing events in 
this system (e.g. Koopmans & de Bruyn 2000; Koopmans et 
al. 2001). We stress the large uncertainties associated with 
these extrapolated calculations. 



9 FUTURE SURVEYS 

As we have discussed, the main reason why it is hard to 
determine the halo density of the observed WD is that the 
sample surrounds the SGP and that we have poor resolution 
of the vertical velocity W . In order to choose between var- 
ious possible local distribution functions, it is necessary to 
repeat the survey at lower Galactic latitude. We discuss two 
potential surveys for the near future that allow a rigorous 
test of some of these underlying assumptions. 



9.1 The Galactic Anti-Center 

In Fig.ll, we show the probability distribution of proper 
motions, seen in the direction of the Galactic anti-center 
(£ — 180° and b = 0°). For simplicity we adopt the same 
survey limits as used by Oppenheimer et al (2001) and a 
standard-candle WD with an average absolute magnitude 
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Figure 11. The probability density of thick-disc plus halo WD 
proper motions towards the Galactic anti-center (see text). The 
contours are spaced by factors of two, decreasing from the maxi- 
mum. The inner dashed circle indicates the proper-motion limit of 
0.33 arcsec yr -1 , whereas the outer dashed circle is indicative of 
3.0 arcsec yr . The extension towards negative fiy exemplifies 
the asymmetric drift. The oblateness of the probability contours 
is related to the ratio oyy /<Tu- 



Mrsqf — 14.3. The latter translates into a maximum sur- 
vey depth of ~125 pc. We used our most-likely model in 
the calculation. If a WD survey in this particular direction 
shows a proper-motion distribution similar to that in Fig.ll, 
it would be strong support to the notion that the WDs with 
high velocities in the U — V plane also have high veloci- 
ties in the VK-direction and therefore form a genuine halo 
population. A significantly flattened proper-motion distribu- 
tion would support a thick-disc interpretation of the sample 
found by Oppenheimer et al. (2001), but requires a rather 
artificial VDF. 



9.2 Halo White Dwarfs with the ACS 

A new opportunity to measure WD proper motions is pre- 
sented by the Advanced Camera for SurveysPI (ACS). It is 
likely that it will be proper-motion limited, rather than mag- 
nitude limited. Scaling from the astrometry carried out by 
Ibata et al. (1999) on the HDF, and requiring that the trans- 
verse velocities be measured to an accuracy ~30 kms" 1 , we 
estimate that a "wide" survey to a depth ~1.2 mag shallower 
than the HDF(N) could be carried out to a limit 



wide 
Mlim 



3.0 dipc x 10 



0.4(M R -15) 



mas yr 



"I http://www.stsci.edu/cgi-bin/acs 
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over Sis ~ 10 4 ster rad (~0.3 square degree) and a "deep" 
survey could reach 



deep 
Mlim 



0.3 4p C xlO°- 4(M «- 15) masyr-\ 



to a depth ~ 1.2 mag deeper than the HDF(N) over S7s ~ 
10~ 6 ster rad (~120 square arcmin). The expected number 
of halo WDs in these fields is 

/f»±/nii: 



dss 2 F h (r,9;v r ,v 4 ,,ve), (25) 



where s measures the distance along the line of sight. Eval- 
uating this integral for the direction of the HDF(N), we find 
that iVwD « 10 ■ (lO°- 3(1B - M « ) ) for the wide survey and 
Nwd ~ 0.6- (10 03(1b - Mr) ) for the deep survey. Evaluating 
the luminosity for the high-velocity WDs, we find that the 
number density per magnitude is roughly constant in the 
range 13<Mjj<16 mag. This gives us 



ji rhalo 
'»WD 



5 (wide) or N v 



0.3 (deep). 



The wide survey will typically observe to a distance of 
~0.5 kpc, whereas the deep survey should see to ~1 kpc. 
The expected number in the existing WFPC2 observations 
of the HDF(N) is ~0.1 and the negative report of Richer 
(2001) is not surprising. 

Of course, these calculations are only illustrative. More 
WDs are expected to be found associated with the thick disc 
and would be seen along lines of sight that look closer to the 
inner Galaxy. However, it is clear that if we can measure the 
white-dwarf colors and proper motions over a few directions, 
then we can use the distribution function approach, outlined 
above, to measure the WD density throughout a major por- 
tion of the halo. 



10 DISCUSSION 

We have analysed the results of Oppenheimer et al. (2001) 
and argued that there appear to be at least two kinemati- 
cally distinct populations of WDs (at the >99% C.L.), which 
one might term a thick disc and a flattened halo population, 
respectively. In contrast, these populations are indistinguish- 
able in their luminosity, color and age distribution. Our 
most likely model indicates that the halo WDs constitute 
~0.8 percent of the nominal local halo density, although this 
fraction decreases as PwD/ptotai~ 0.008 • (r/r c ) _1 '° with in- 
creasing Galacto-centric radius r/r c , for an assumed Galac- 
tic density profile ph oc r~ 2 . A more conservative lower limit 
of 0.3 percent (90% C.L.) can be placed on the local halo 
mass fraction in WDs. This WD density is 5.0tfl (90% 
C.L.) times higher than previously thought and comparable 
to the local stellar halo density (Gould et al. 1998). This 
is in conflict with the estimates of WD densities, based on 
population synthesis models and Salpeter initial mass func- 
tions, and requires a non-standard explanation if one re- 
quires these WDs to be formed in the halo. Based on the 
notion that 90% of the WDs, found by Oppenheimer et al. 
(2001), arc proper-motion limited (Sect. 6), we do not expect 
the white-dwarf luminosity function to rise sharply beyond 
the survey magnitude limit (i.e. WDs with ages >10 Gyr), 
although a distinct very faint population of WDs with ages 
3>10 Gyr, which are not related to the presently observed 
population, can not be excluded. 



We find that the halo WD population is flattened (q = 
(c/a)p ~ 0.9), and expect q to become somewhat smaller 
once a more realistic flattened potential of the disk, bulge 
and halo is used in the models. The value for q that we quote 
should be regarded as an upper limit. 

The microlensing optical depth inferred from the halo 
WD population is estimated at t^d ~ 1.3 x 10 -9 . Even if 
we include the contribution from the thick disk, the opti- 
cal depth remains Z-,5% of that observed toward the LMC 
(tlmc ~ 1-2 x 10~ 7 ; Alcock et al. 2000). We conclude that 
the population of halo WDs, observed by Oppenheimer et al. 
(2001) and on orbits that are probed by the velocity distri- 
bution function in the Solar neighborhood, do not contribute 
significantly to the observed microlensing events toward the 
LMC. One avoids this, if either ~25 times more WDs are on 
orbits with perigalacticons well outside the Solar neighbor- 
hood (for example a mass shell) or on resonant orbits. Both 
types of orbits could be undersampled locally. Similarly, as 
mentioned above, if there exists yet another population of 
very faint old (~15 Gyr) WDs in the halo, they could also in- 
crease the number of microlensing events towards the LMC. 
However, we regard both solutions as unlikely and artificial. 
In addition, both solutions would violate constraints from 
metal abundances (see Sect.l) and (simple) population syn- 
thesis models. In Sect. 8, we showed that a similar "shroud" 
of white dwarfs around the LMC itself, in principle, can ac- 
count for the observed optical depth, although the result is 
strongly model dependent. 

We have performed a number of consistency checks of 
the data. For example, we find excellent agreement between 
the maximum-likelihood value of the asymmetric drift, the 
statistical value and the value inferred from other stellar 
population in the solar neighborhood. There is also no ev- 
idence for either a drift of the WD population in the U 
or W direction, which could have indicated that the pop- 
ulation was not in dynamic equilibrium or could be part 
of a stellar stream or group. In both cases our approach, 
based on the collisionless Boltzmann (or Vlasov) equation, 
would have been invalid. In addition, the thick-disc likeli- 
hood results agree well with the empirical relation between 
the asymmetric drift (v a ) and the radial velocity dispersion 
((Tu), determined from a range of local stellar populations. 
We also find that the asymmetric drift of the halo popula- 
tion is close to that (i.e. v a ~ v c ) expected for a population 
that is mostly pressure supported (i.e. au ~ v c /\/2). If we 
assume that o\j=oy for the WD population, the likelihood 
decreases by a factor ~ 10 6 and the thick-disc result no 
longer agree with the empiral relation. A flattened velocity 
ellipsoid (au<crv) naturally explains the relative scarcity of 
WDs in the regions around (U, V) = (0, —220) kms -1 and 
V < -220 kms -1 , noticed by Reid et al. (2001). 

The absence of any clear correlations of velocity with 
color, absolute magnitude or age, and the dynamical self- 
consistency checks give us confidence that the population of 
WDs found by Oppenheimer et al. (2001) does not contain 
a significantly hidden bias. Especially, the fact that the sur- 
vey is 90% proper-motion limited, not magnitude-limited, 
underlines that we do not miss many WDs that are too 
faint (but see the comments above). The range of luminosi- 
ties should thus provide a reasonable representation of the 
underlying luminosity function of this particular population 
of halo WDs. 
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These internal and external consistency checks should 
make our results quite robust. However, there are uncer- 
tainties associated with errors in the photometric parallaxes 
which feed into both the densities and velocities. In addi- 
tion, several of the WDs in this sample might belong to 
binaries, although an inspection of the two-color diagram 
(Fig.l) suggests that also this is unlikely to be a major prob- 
lem. Only a systematic overestimate of distances by a factor 
<;2.5 would move most high-velocity WDs into the thick-disc 
regime, but would also result in a typical velocity dispersion 
of the sample only half that of the thick disc and average 
WD luminosities that are too small. The absence of appar- 
ent correlations with the WD velocities, makes it unlikely 
that distance overestimates would preferentially occur for 
the highest velocity WDs. Another way out of the conclu- 
sion that there exists a significant population of halo WDs, 
would be to accept that the velocity ellipsoid of these high 
proper-motion WDs is significantly flattened in the vertical 
direction. Although, this would be consistent with the obser- 
vations, as we indicated in Sect. 5, it would require a highly 
unusual VDF, unlike what is observed for both young and 
old stellar populations in the Solar neighborhood. 

We have also provided further evidence that the fast 
and slow WDs come from the same population. This would 
not be a surprise if they were mostly born in the disk and 
then deflected dynamically to high altitude. We have pro- 
posed a possible mechanism that could preferentially eject 
WDs from the disc, involving orbital instabilities in evolving 
multiple (e.g. triple) stellar systems (§7.2). We showed that 
the total halo plus disc WD mass of the Galaxy (i.e. ~15% 
of its stellar mass) is roughly consistent with that expected 
from the total stellar mass in our Galaxy, based on stan- 
dard population synthesis models. The agreement suggests 
that not a re-thinking of Galactic starformation models is 
required, but that the real question is how to eject disc WDs 
into the halo with high velocities. We propose that this can 
be achieved through the orbital instabilities in evolving mul- 
tiple stellar systems. 

The key to understanding the WD distribution with 
more confidence is undoubtedly to measure the component 
of velocity perpendicular to the Galactic disk and this is best 
accomplished with a proper motion survey at lower latitude. 
We have suggested two possible surveys for the near future: 
one towards to Galactic anti-center, similar to the survey 
discussed in this paper, and one that could be done with 
the Advanced Camera for Surveys (ACS) on HST. The first 
survey could unambiguously show that the velocity disper- 
sion perpendicular to the disc is comparable to that parallel 
to the disc. Measuring the density of WDs as a function of 
Galacto-centric radius (for example with the ACS) , together 
with a local WD density calibration, can provide strong con- 
straints on the shape of the Galactic potential, as is apparent 
from Eqns (17) and (18) in Sect. 7. These constraints im- 
prove dramatically if we can also constrain the VDF at the 
same points in the halo. This is especially relevant in case 
of the halo WDs, even though they are faint, because their 
high velocity dispersion allows one to probe large (tens of 
kpc) Galacto-centric distances and still be in dynamic con- 
tact with the Solar neighborhood. The relation between the 
VDFs at different radii is then given through the integrals of 
motion that are conserved along orbits between these points 
and the Galactic potential. Measuring the density and VDFs 



along these orbits can then, in principle, be used to recon- 
struct the Galactic potential. It is also surprising that the 
scale height of WDs, kinematically belonging to the thin 
disc, appears at least twice that previously thought (Ma- 
jewski & Siegel 2001). Could this population form a bridge 
to a flattened white-dwarf halo population? 

We conclude that the discovery of a surprisingly large 
population of high velocity, old WDs, made possible by ad- 
vances in understanding their atmospheres and colors, is a 
significant accomplishment and seems to opens up a new un- 
expected window into the stellar archaeology of our Galaxy. 
However, to confirm or reject the ideas put forward by Op- 
penheimer et al. (2001) and in this paper a significantly 
large sample is required, not only out of the Galactic plance, 
but also in the Galactic plane. In addition deeper surveys 
(e.g. with the ACS) could probe much farther into the halo 
and possibly even detect WDs that are significantly fainter. 
Deeper surveys would also be able to probe the transition 
between thick-disc and halo and be less 'contaminted' by 
disc WDs. 
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